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P\| ' Abstract 

One considers a class of PDEs systems (1) equivalent to (8) and one 
►-^ , determines the PDEs system (10) which defines the associated symmetry 

fvq ' group. Particularly, for the Blair system (2)+(3) equivalent to (14), one 

^H ' finds the symmetry Lie group G (Theorem 5). One proves that the class 

C^ , of PDEs systems in the form (8) + (14') which are invariant with respect 

f^ ' to the Lie group G, is reduced to the Blair system (Theorem 6). One finds 

new solutions B\ and BJ of the Blair system, and thus new "force- free" 
model of solar physics. 






X 



Key-words: Blair system, infinitesimal symmetries, symmetry group, group- 
invariant solutions. 

Mathematical Subject Classification: 58G35, 35A30 
Physics Subject Classification: 96.60.H 

1 Introduction 



The classical method for finding the symmetries of PDEs reduces to the Lie 
5^ ' group method of infinitesimal transformations. The classical symmetries of 

5y \ PDEs are related to vector fields which are functions of the independent and 

dependent variables. A symmetry group is a Lie group which transforms the 
solutions of the system into itself, thus it provides a mean of classifying dif- 
ferent symmetry classes of solutions, where two solutions are considered to be 
equivalent if one can be transformed into the other by the same of the group 
element. The infinitesimal generator is a vector field on the underlying mani- 
fold which determines a fiow (1-parameter group of transformations). One can 
regard the entire group of symmetries as being generated in this manner by the 
composition of the basic fiows of its infinitesimal generators. In [6] G.Bluman 
and J.D.Cole introduced the algorithm of finding the solutions of the system 
by using the transformations of symmetry. A modern presentation of this the- 
ory can be found for example in [12], [15], [19]. There are many applications 



of this theory in the study of the PDEs systems which arises of mathematics, 
mechanics and physics, for example [l]-[4], [6]- [19]. 
The solutions of the vector equations 

(1) curl B = f -B 
and 

(2) div B ^0, 
where / is an arbitrary function, and 

d d d 

B = u{x,y, z)-- + v{x,y,z)—- +w{x,y,z)--, 
ox oy az 

is a differentiable vector field on a simply connected domain D C R^, define 
"force-free" model of solar physics. On the other hand, the vector equation 

(3) curlB = |B|-B, 

was introduced by Blair in [5] : a solution of it gives a conformally flat contact 

metric structure on R^. The vector fields which satisfy the vector equation (3) 

are 

r, ■ 9 d 

i>i = smz— h cosz— — , 

ox oy 

and 

8(xz-y) d 8(x + yz) d 4{1 + z'^ - x'^ - y'^) d 



B- 



(1 + x^ + y^ + z^Y dx (1 + a;2 + y2 _|_ ^2~)2 Qy (^l + x'^ + y^ + z^)^ dz ' 



and Bi is a solenoidal vector field. In the cylindrical coordinates, by using the 
method of succesive aproximations, Blair gives a solution of the equation (3) 
for which (2) is also satisfied, thus a new solution for the "force- free" model 
equations. The finding of the PDEs system (2)+(3) solutions is still an open 
problem. 

In this paper we shall determine infinitesimal symmetries associated to the 
PDEs system (1) for / = f(u,v,w) (Theorem 3) and in the particular case for 
f ^\B\ (Theorem 4). We cah Blair system the PDEs system (2) -I- (3). Also we 
shall find the symmetry group G associated to it (Theorem 5). We shall prove 
that the only PDEs system (l)+(2) which is invariant with respect to the group 
G is the Blair system (Theorem 6). We shall prove that the known solutions 
are group-invariant solutions and by using the symmetry group G one finds new 
solutions of it: B[ ' and B[ . 

We make the remarks that there are many computational symbolic programs 
for finding the defining system of infinitesimal symmetries and ours calculus was 
verified by using the Head's program LIE [10]. This paper gives a new point of 



view for the PDEs systems which appear in the solar physics and as well in the 
contact geometry. 

We adopt the notations of the Olvcr's book [12]. We start to make a short 
presentation of the symmetry group theory in the general case of the PDEs 
system. 

2 Symmetry group for PDEs system 

Let us consider the PDEs system 

(4) A,(x,u(")) = 0, iy = l,...J, 

with X = (x^, ..., a;P) the independent variables, u — [u^, ...,u'') the dependent 
variables and 

A(x,k(")) = (Ai(a;, «(")),. ..,Ai(x,w("))), 

where we denote u*-"-* all the partial derivatives of the function u to up n. Any 
function u = h(x), 

h-.DcW^UcR.'^, h^{h\...,h''), 

induces the function 

called the n-th prolongation of h, which for each x E D, pr^^^>h{x) is a vector 
whose qp^'^' = qCp^„ entries represent the values of h and all its derivatives up 
to order n at the point x. 

The total space D x U^'^\ whose coordinates represent the independent 
variables, the dependent variables and the derivatives of the dependent variables 
up to order n, is called the n-th order jet space of the underlying space D x U. 

Thus A is a map from the jet space D x t/^"^ to R'. The PDEs system (4) 
determines a subvariety 

of the total jet space D x U'-"'K One identifies the PDEs system (4) with its 
corresponding subvariety S. 

Let M C D X U he an open set. We suppose X is a vector field on M, with 
corresponding local 1-parameter group exp(eX). The n-th prolongation of X, 
denoted by pr'^^-^X, is a vector field on the n-jet space M^") , and is defined to be 
the infinitesimal generator of the corresponding prolonged 1-parameter group 
pr(")[exp(£X)]: 



d 

(x,m(")) — 

for any (a;,u(")) G M(") 



Pr-(")X|(,^„(„,) = -pr(")[exp(£X)](x,«("))U=o 



The PDEs system (4) is called to be of maximal rank if the Jacobian matrix 



JAix,u^">) 



(«)^ - 



of A, with respect to all the variables (x, u*^"^), is of rank I whenever 
Theorem 1. Let 

i—l a—1 

be a vector field on open set M <Z D x U . The n-th prolongation of X is the 
vector field 

(5) p^(")x = X + ^5]a>;^(x,z.("))-^, 

defined on the corresponding jet space M^"' (Z D x U^'"'' , the second summation 
being over all multi-indices J — {ji,---jk) with 1 < jk l£ P, ^ l£ k < n. The 
coefficient functions $;^ ofpr^"'X are given by the following formula 

$f (x, u(-^) ^ Dj U^ - J2 C<] + E CW},' 

\ i=l / i=l 

where uf = |4, <, = #4- 

Theorem 2. {Infinitesimal criterion of invariance): Let us consider the 
PDEs system (4) of maximal rank defined over M d D x U . If G is a local 
group of transformations acting on M , and 

(6) pr(")X[A^(x,u("))] =0, v = l,...,l, 

whenever A,j{x,u^''^') = 0, for every infinitesimal generator X of G, then G is 
a symmetry group of the system. 

Proposition 1. If the PDEs system (4) defined on M C DxU is of maximal 
rank, then the set of infinitesimal symmetries of the PDEs system forms a Lie 
algebra on M. Moreover, if this algebra is finite-dimensional, then the symmetry 
group of PDEs system is a Lie group of local transformations on M . 

Algorithm for determination of the symmetry group G associated 
to the PDEs system (4): 

-one considers the field X on M and one writes the infinitesimal invariance 
condition (6); 



-one eliminates any dependence between partial derivatives of the functions 
u", determined by the PDEs system (4); 

-one writes the condition (6) like polynomials in the partial derivatives of 
u"; 

-one equates with zero the coefficients of partial derivatives of u" in (6), 
written as polynomials in the derivatives of the functions u"; it follows a PDEs 
system with respect to the unknown functions C*, 0" and this system defines 
the Lie symmetry group G of the given PDEs system. 

In general, for each s-paramctcr subgroup H of the full symmetry group G 
associated to the PDEs system {in p > s independent variables ), there will 
correspond a family of group-invariant solutions. Thus, a classification of these 
solutions is obtained by using an optimal system of group-invariant solutions 
from which any other solution can be derived. 

Proposition 2. If u — h{x) is an H -invariant solution of the PDEs system 
and g G G, then the composed function v = h{x) = g ■ h{x) is a H -invariant 
solution, where H ~ gHg^^ is the conjugate subgroup under g. 

The problem of classifying group-invariant solutions reduces to the problem 
of classifying subgroups of the full symmetry group G under conjugation and 
this is equivalent with the classifying subalgebras of the Lie algebra g of the 
group G under the adjoint representation. 

An optimal system of s-parameter subgroups is a list of conjugacy inequiva- 
lent s-parameter subgroups with the property that any other subgroup is conju- 
gate to precisely one subgroup in the list. Similarly, a list of s-parameter subal- 
gebras forms an optimal system of s-parameter subalgebras if every s-parameter 
subalgebra of g is equivalent to a unique member of the list under some element 
of the adjoint representation. 

Thus we compute the adjoint representation Ad G of the underlying Lie 
group G, by using the Lie series : 

°c „ 2 

(7) AdiexpieX)Y) = ^ _(adX)"(r) = Y-e[X,Y] + -[X, [X,Y]] - ... 

ri=0 

3 Symmetries of PDEs systems in solar physics 
and contact geometry 

Let us consider 

d d d 

B ^u{x,y,z)— + v{x,y,z)— +w{x,y,z) — , 

a differentiable vector field on a simply connected domain D C R'^. Let 
(8) 



Wy~Vz - 


= uf 


Uz — Wx = 


- vf 


Vx~Uy -- 


- ^f- 



be the PDEs system associated with the vector equation (1), where / is an 
arbitrary difFerentiable function of u, v, w. We denote 

Ai ^Wy -v^- uf, A2 ^ u^ - Wx - vf, A3 ^ Vx - Uy - wf, 

and compute the partial derivatives 

au aw ow ov^ owy 

aA2 ^ 9A2 ^ ^ aA2 ^ aA2 _, aAa 
^— = -f /«, ^— = -J - w^, ^— = -«/«,, ^ — = 1, ^ — = -1, 

ou ov ow ouz owx 

dAa , 9A3 , aAg _ _ 9A3 ^ 9A3 ^ 

OU ov ow ouy ovx 

Let Ja be the Jacobi matrix of the function A. It results that rankJ^ ~ 3 
and thus the PDEs system (8) is of maximal rank. Let us consider 

^ d d ^ d ,9 ,9 d 

ox oy Oz ou ov ow 

a vector field on an open set M C D x U oi the space of the independent 
and dependent variables of the system, where (^,ri,9, (f>, X^i}) are functions of 
x,y, z,u,v,w. 

If we consider X to be the infinitesimal generator of the symmetry group of 
PDEs system (8), then the first prolongation of it is the next vector field 

(9) pr^^^X = X + $"-^ + $^-^ + $"^ + A"^ + A^'^ + A^^+ 

OUx OUy OUz OVx OVy OVz 

d d d 

+^^^- — + my- — + m- 



dwx dwy dwz ' 



where 



$^ = 4>x + Ux{4>u - Cx) - Uy-qx - Uz9x + <i)vVx + (t'wWx - ulQu- 
Ux^y^u ^x^z^u '^x'^xsv ^y^x^v '^x^z^v ^x^x^w 
- UyWxT]^ -UzWxO^, 

$2/ ^ (jjy - UxCy + Uy{(j)u - Vy) - Uz6y + Vy4>y + Wy4>,u] - UxUyCa - UyTju- 

^- liylLzyJu UxVyi^y — UyVyTjy ^^ Uz^y'^V ^^ ^^x'^y^W ~~ UyWyTJw ^~ UzWyUw-) 

$^ = (t)z -UxCz-Uyr]z+Uz{(f)u -Oz) +Vz(f)v+Wz(f)w -UxUzCu-UyUztJu- 

- ujOu - UxVzCv - UyVzrjy - UzVzOv - UxWzCw - UyWzTjw - UzWzOw, 
A^ = K + UxK + VxiK - Cx) - VyVx - VzOx + WxXw - UxVxCu - UxVyTlu- 

- UxVzOu - vlCv - VxVyT]y - VxVzOv ^ VxWxCw ~ VyWxTjw -~ VzWxO^, 

6 






UyVyTju - UyVzOu - V^VyQv - VyTj^ - VyVz9y - V^WyCu 



- VyWyT]^ -V^Wytl^, 

A^ = X^+U^Xu-V^Cz-VyTI^+V^X^v -0z) +W2XW -VxU:,(u-VyU^r]u- 

- VzU^Ou - V^W^Cw - VyW^T]^ - V^W^e^ - V^V^Cv - VyV^Tj^ - V^O^, 
*^ = Ipx + U^tpu + V^^py + W^{tp^ - C^) - Wyl]^ - W^9^ - U^W^Cu- 

- u^Wyiriu - u^w^9u - v^w^Cy - v^Wyrjy - v^w^O^- 

- wlC^yj - W^WyT]^ - W^W^Oy,, 

"^y = IPy +Uy'4)u+Vy1py -WxCy +Wy{lpn] -Vy) -Wz9y -WxUyQ- 

- WyUyTJu - W^Uy9u - W^VyCy - WyVyTly - W^Vy9y " W^Wy^^- 

- wlr]y,~WyW^0y,, 

*^ = 1pz+UzTpu+Vz1pv+Wz{lpw- Oz) -WxCz-VzWy -UzWxCu-UzWyr]u- 

- 6uUzWz - VzWxCv - VzWyTJy - VzWzOy " W^W^Cw " WzWyTJyj - wl9yy. 

In this case, the infinitesimal invariance condition (6) impHes 

-{f + uU)^-ufyX-uf^^i, + -^y-K' = 
-vfu^-{f + vfy)X-vU^ + <^''-^^ = 

~wM~wfyX-{,f + wU)ij + K^-^y = 0. 

Substituing the functions $^, (f>^, A^, A^, ^^, ^i* and after ehminating any 
dependencies among the derivatives of the u,v,w caused by the PDEs system 
(8) itself, we find the following PDEs system which defines the symmetry group 
of the studied PDEs system 

■4^y-'Xz~'f>{ufu+f)-uXfy-uipJy,-vfXu+wfC,z+uf{il;y,~riy)+vwpCu-u^.f'^riy,+ 

+ Uy{ll;y+Q+VfCy-Ufr]u)+Vy{lPy+f]z+VfTJy-UfT]y)+Vz{9z-Xy+tP^-TJy + 
+ Vf9u+WfCy-2uJriyy)-Wx{Xu+C,y-WfCu+UfC,y,)-Wz{Xyy+9y-WfCy,+Uj9yy) + 

+UyWz{Cw-0u)+VyWz{Vw-0y)+UyVz{Cv-Vu)+WxVz{9u-Cw)+VyWx{r]u-Cy) =0. 

4>z-'4'x-v fu4>-{v fy+ f)X-v fw'4}+v f {4>u-9 z)-w fijjy+ufrjx-v^ f 9u+UW frfy- 
-Ux{Cz+'(pu+v.fCu-ufriu)-Uy{riz+'ipy+vfr]u-ufT]y)+Vz{(j)y+r]x-vfdy+wfr]y)+ 
+Wx{Cx -■>pw + 4'u-9z- 2vf6u + wfC,y + ufrjw) + Wz{(t>w + 9x- vfOy, + Wf9y)+ 

+UxVzir]u - Cv) + Ua:Wz{9u - (w) + UyWz{9y - r]y,) + wliC^yy - du) + WxUy{Cy- 

-rju) + WxVzi-qw - 9y) = 0. 

K~4>y~W(t)fu-WXfy-iwfy,+f)lP+Vf9y+Wf{Xy-Q)-Uf4)y,-w'^PCy+UVf'^9^ + 



+Ux{Xu + Cy + ufCw - wfCu) + Uy{riy - (t>u + At, - Cx + vfdu + ufrjw - 2wfCv)- 
--Vy{(f)y + r/a; - vf9^ + w f Tj^) - V ^ {9 ^ +(/>«,+ wf9^ - vfd^) + Wx{X^ + 9y- 

-wfCni+uf9^)+ul{T]u - Cv)+Ua;Vy{Cv - Vu) + U^V^iCw - 9u) + WxVy{9y - T]^) + 

+ WxUy{9u - Cw)+ UyVziVw ^ 6 v) = 0- 

If we equate the coefficients of the remaining unconstrained partial deriva- 
tives of u,v,w to zero, we obtain the following PDEs system which defines the 
symmetry group of the studied PDEs system: 

V'y - ^2 - (/>(u/„ + /) - uXf^ - U^Jw - vfXu + wfC,:,+ 

+uf{ipw -Vv) + vw,fCu - W^.fllw = 

(f>z-'(px- vfu4' - (vfv + /)A - vf^ip + Vf{(f)u - 9^)- 

-wftpy + ufrjx - v^f^9u + uwf^rjy = 

\X - (f>y~ W(f)fu - WXfy - {Wfy, + f)lp + Vfdy + W f {Xy - (x) ~ 
-Uf(t>y, - w'^f^Cy + UVp9y, == 

01) = -rjx + fivVw - WT]y) 
(10) 0,. - -9x + f{v9^ - w9y) 

K = -Cy + .f{wC,u - U(yj) 

A«, = -9y + f{w9u - u9w) 
i^u = -Cz + .f{uCy - v(u) 
ipv = -Vz + fiuriv - vrju) 

(l^u-lpw = -Cx +0z+ fC^vOu - wCd - UT]yj) 

i'w ~ K = -dz +% + fC^Uriy, - Vdy - WC,y) 

Vu — Sv 

^w — "li- 
lt results the next theorem 

Theorem 3. The general vector field which describes the algebra of infinites- 
imal symmetries associated to PDEs system (8) is 

B d „ d d d d 

ox oy oz ou ov aw 

where the functions Q, rj, 6, (j), X and -0 satisfy the PDEs system (10). 



Wy - 


- Vz 


— U^V? + u^ + w^ 


Uz - 


-Wx 


= Vyju?- + U^ + u;2 


Vx - 


- Uy 


— w\/u^ + w^ + w^. 



Now let us consider the PDEs system 

(11) 

associated to the vector equation (3). 

We obtain 

Theorem 4. A Lie algebra of infinitesimal symmetries associated to the 
PDEs system (11) is described by the next vector fields 



(12) 



^1 = 




d 

dx 


d 
dy 


d 

du 


d 

dv 


X2-- 




d 

dy 


■ + y 


d 
dz 


— u 


dv 


d 

; 

dvu 


^3 


= - 


d 

ox 


d 

Oz 


d 

- w- [ 

du 


d 

dw 


X, 


_ = 


d 

dx' 


X5 


= 


d 

dy' 


Xe^ 


d 
dz' 






Xr 


d 


d 

dy 


d 

dz 


■ — u 


d 

du 


— V 


d 
dv 


- w 


d 

dw' 








Xs 


; = 2XZ' 


d d 

dx dy 


{-"- 


x^- 


v'\ 


'1 


+ 2(: 


xw 


- zu}- — h 
du 





d d 

+2{yuj — zv)— 2(xu + yv + zw)— — , 

dv ' dw 

9 ^ , 2 2 2\ 9 d . . d 

X9 = xy—- + -(y -X - z )—- +yz—- + +{xv- yu)- 

dx 2 dy dz du 

I \ 9 , , d 

— [xu + yv + zw)- — h (zv — yw)——, 
dv dw 

1 d d d d 

-'^10 = 77(2;^ - y^ - ^^)Tr + ^y^r + ^^^ — (^" + yv + zw)- — h 

2 dx dy dz du 

, \ 9 , . d 

+ [yu — XV)— — h (zu — xw)— — . 
dv dw 



Proof. We substitute / — \Ju^ + w^ + w^ in the PDEs system (10). If 
we consider ^, ry, 0, </), A and ^/^ polynomial functions of x, y, z, u, u, w, we get the 
following solution: 

( = C7X -Ciy-C3Z + 2(78X2 + Cgxy + iCio(a;2 -y^- z^) + C4 
j^ ^ Cix + C^y - C2Z + 2Csyz + ^Cg{y^ - x'^ - z^) + Cioxy + C5 
9 = C3X + C2y + +C\z + Cs{z'^-x'^-y'^) + Cgyz + Cioxz + Ce 
(j) = —Cju — Civ — C3W + 2C%{xw — zu) + C(){xv — yu) — 
(13) —Cii^{xu-\-yv + zw) 

A = Ciu — C7U — C2W + 2Cs{yw — zv) — Cq{xu + yv + zw)+ 

+Cio{yu — xv) 
ip — C3U + C2V — C7W — 2Cs{xu + yv + zw) + C9{zv — yw)+ 

+Cio{zu — xw), 



with Ci e R, i = 1, ..., 10. The infinitesimal generator of the associated sym- 
metry subgroup is X = J2i=i CiXi, where the vector fields Xi are given by the 
relation (12). 

The structure constants of the Lie algebra generated by the vector fields (12) 
are given by the next table 



[■,•] 


Xi 


X2 


X3 


Xi 


Xs 


Xe 


X7 


Xs 


X9 


Xio 


Xi 





X3 


-X2 


-Xs 


Xi 











Xio 


-X9 


X2 


-X3 





Xi 





-Xe 


Xs 





2X9 


~kXs 





Xs 


X2 


-Xi 





-Xa 





Xi 





2X10 





-i^s 


Xi 


X5 





Xe 











Xi 


-2X3 


-Xi 


X7 


Xs 


-Xi 


Xe 














Xs 


-2X2 


X9 


Xi 


Xe 





-Xs 


-Xi 











Xe 


2X7 


X2 


X3 


Xr 











-Xi 


-Xs 


-Xe 





Xg 


X9 


Xio 


Xs 





-2Jf9 


-2X10 


2X3 


2X2 


-2Xy 


-Xs 











Xg 


-Xio 








Xi 


-Xg 


-X2 


-Xg 











Xia 


Xg 





kXs 


-Xj 


-Xi 


-Xs 


-Xu, 












For the Blair system 



(14) 



Ux +Vy + Wz 



= uf 
= vf 
= wf 



where 
(14') 



Ux+Vy+Wz= 0, 



is equivalent to (2), it results 

Theorem 5. The following vector fields 



^^ d d d d 

X\ = -y^r- + X- V—- + u—-, 

ox oy ou av 



(15) 



d d d d 

oy oz Ov ow 



X. = -z 



d_ 
dx 



d 



d 



X. 



dx' 



X. 



d 



X- w—- +u—-, 

oz ou ow 

— , Xe^—. 
dy' dz 



10 



Xr 



dx 



d_ d_ 

dy dz 



d 



d 



du dv 



w 



dw' 



describe the Lie algebra g of the infinitesimal symmetries associated to the Blair 
system (14)- 

Proof. The PDEs system (14) is of maximal rank. By using the above 
algorithm, wc find that 



(16) 



where Ci G R, i = 1, •.., 7, is the solution of the PDEs system which defines 
the symmetry group of Blair system. It results that X = X]i=i CiXi, is the 
infinitesimal generator of it, where the vector fiels Xi are given by (15). 

Remark. The Lie algebra g associated to the symmetry group of the Blair 
system is a subalgebra of the algebra of infinitesimal symmetries associated to 
PDEs system (11). Using the Lie series, one determines the adjoint representa- 
tion AdG of the Lie symmetry group G: 



c = 


= Crx - Ciy - C3Z + d 


V = 


= C\x + Cjy - C2Z + C5 


e -- 


= Czx + Cay + +C7Z + Cg 


<P -- 


= —Cju — Civ — C3W 


A . 


- C\u — C-jv — C2W 


^ = 


= Czu + C2V - C7W, 



Ad 


Xi 


X2 


X3 


Xi 


Xi 


X2 cos £ — X3 sin e 


X3 cos e+X2 sin e 


X2 


Xi COS S+X3 sin 6 


X2 


X3 cos e — Xi sing 


X3 


X-i cos S — X2 sin s 


X2 cos s+Xi sin s 


X3 


Xi 


Xi-eXs 


X2 


X3-eXe 


X5 


X^+eXi 


X2-eXe 


X3 


Xe 


Xi 


X2+eXr, 


X3+eXi 


X7 


Xi 


X2 


X3 



Ad 


Xi 


X5 


Xe 


Xr 


Xi 


Xi cos s+X^ sin e 


X5 COS €—Xi sin € 


Xe 


X7 


X2 


Xi 


Jf 5 cos e+Xe sin e 


Xe cose-Xs sine 


Xr 


X3 


Xi COS e+Xe sin e 


Xs 


Xe cos e — Xi sin e 


X7 


Xi 


Xi 


Xs 


Xe 


Xr-eXi 


Xs 


Xi 


X5 


Xe 


Xr-eXs 


Xe 


Xi 


X5 


Xe 


Xr-eXe 


Xr 


e'Xi 


e^Xs 


e^Xe 


X7 



11 



Remark. If {u = f{x,y,z), v = g{x,y,z), w — h{x,y,z)) is a solution of 
the Blair system, then the following triples 

u^^' = af{ax + by,—bx + ay,z) — bg(ax + by,—bx + ay,z) 
v^^' — bf{ax + by, —bx + ay, z) + ag(ax + by, —bx + ay, z) 
w^^' = h{ax + by,—bx + ay,z), 

u^^i = f{x,ay + bz,—by + az) 

v^"^' — ag{x,ay + bz,—by + az) — bh{x,ay + bz, —by + az) 

w^^i — bg{x,ay + bz,—by + az) + ah{x,ay + bz,—by + az), 

u^^> ~ af{ax + bz,y,—bx + az) — bh{ax + bz,y,—bx + az) 

v^^' = g{ax + bz,y,—bx + az) 

w^ = bf{ax + bz,y,—bx + az)+ah{ax + bz,y,—bx + az), 

{u(4) ^ f{x-e,y,z) ( m(5) == f{x,y-e,z) 

v'-^^ = g{x-£,y,z) I «('"') = g{x,y-e,z) 

w*-^-* = h{x — e,y,z), [ w^^' — h{x,y~e,z), 

u^^' — f{x,y,z~e) ( u"' — e^^f{e^^x,e^^y,e~'^z) 

^ — .9(2^7 2/; -2 ^ s) \ "^ — e~^g{e~'^x,e~^y,e^^z) 

w^^' — h{x,y,z — e), [ w"' — e''^h{e~'^x,e~^y,e^^z), 

are also solutions of the Blair system, where £ G R, cos e — a and sin e ~ b. 

By using the adjoint representation AdG of the grup Lie G we can find an 
optimal system of s-parameter subalgebras associated to the Blair system, with 
s = 1, 2. For example, it results six 2-dimensional optimal subalgebras, which 
are generated by the next vector fields: 

(18) ^4,^5; Xz,Xq; X4,Xq, 
and respectively 

(19) Xi,X(i; X2,Xr,; X3,X4. 

1. Let us consider the symmetry subgroup for which the infinitesimal gen- 
erators are the vector fields X4 and X5 and 

F — G{z,u, V, w) 

the invariant function of it. By substituting 

u — g{z), V = h{z), w — k{z), 

in the Blair system, we get the following system of differential equations: 

fc = 
-h' = g^/ffTh? 
g' = /iVT + T?. 
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By using a coordinate transformation, the solution of it is 

u{z) — sinz, v{z) — cosz, w — 0, 

and it results the vector field Bi. 

2. Now we determine the grup- invariant solutions for the symmetry subgroup 
with the infinitesimal generators Xi §i Xq. The invariant function is 



F = G( v x^ + j/^, xu + yv, XV — yu, w). 
If we consider 



xu + yv = giyx"^ + y"^), xv — yu ^ hi^\/ x'^ + y^), w = 'y{\/x'^ + y"^), 
the Blair system turns in 

.9 = 



^7' = W(^) +7^ 



h' = nJ{^y+i^ 



where r = \/ x"^ + y^, or in an equivalent form 



-i = ii^f W^ 
i/3 + /?' = 77^^+7^, 

by using the change of function /3 = -. The solution of this system was deter- 
mined by Blair using the method of succesive aproximations [5] . 

By circular permutation of x, y, z, one determines other solutions of the Blair 
system. 

According with the above remark, using the relations (17), one finds new 
solutions of the system. For example in the case of the vector field 

r, ■ d 9 

Bi = smz— — h cosz— — , 
ox oy 

the following vector fields are also solutions 

Sf ) = sin(z -e)-^+ cos(z - e)-^ = sf \ 
ox oy 

o o o 

B\ — sin(az — by)— h acos(az — by)— — h fecosfaz — by)-^^ 

ox ' oy oz 

o o o 

B\ — asinfaz — bx)— — h cosfaz — bx)— — h 6cos(az — 6x)— -, 
ox oy oz 
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B\ ' = e~^ sin(e^'^z)— h e^^ cos(e~^z)— — . 

ox ay 

Because B\ ' = B{' = Bi, we get the next result: the solution Bi is invariant 
with respect to the 2-parameter subalgebra described by X4 and X5. Analo- 
gously, one find new solutions of he Blair system by using the Blair solution. 
We make the remark that the vector field 

8(xz-y) d 8(x + yz) d 4{1 + z'^ - x'^ - y'^) d 

B-2 — ' 



(1 + a;^ + J/2 + z'^Y 9x (1 + x^ + y^ + z^)^ dy {1 + x^ + y^ + z^y dz ' 

which is a solution of the vector equation (3), can be found by using the subgroup 

generated by Xi, because B2 = -62- 

Now we study the inverse of the Theorem 5. 

Theorem 6. The only PDEs system (8) + (14' ) which is invariant with 
respect to the symmetry group G associated to the Blair system is the Blair 
system. 

Proof. One writes the infinitesimal conditions (6) in the case of the PDEs 
system (8)-|-(14') and one substitutes the vector fields given by the relation (15). 
It results the next PDEs system 

ufu + vfv + wf^ = f 

ufv - vfu = 

vfw - wfy = 

ufw - wfu = 0, 



with the solution f{u,v^w) = \/ v? + f^ + w^. 

We consider that the study of the PDEs system (l)+(2) from the point of 
view of the symmetry group theory is very interesting, because for this PDEs 
system are many known solutions, but we don't have any classification of them. 
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